Dynamics of Vortex Shells in Mesoscopic Superconducting Corbino Disks 



o 
o 

(N 

■4— » 

o 
O 

o 

m 

"a 

o 
o 

Oh' 
o3 



13 

o 
o 



> 

m 
oo 

o 

O 

-i— > 
c3 



i 

c 

o 
o 



X 



V.R. Misko and F.M. Peeters 
Department of Physics, University of Antwerpen, 
Groenenborgerlaan 171, B-2020 Antwerpen, Belgium 
(Dated: February 6, 2008) 

In mesoscopic superconducting disks vortices form shell structures as recently observed in Nb 
disks. We study the dynamics of such vortices, driven by an external current Io, in a Corbino 
setup. At very low 7o, the system exhibits rigid body rotation while at some critical current 7 Cj i 
vortex shells rotate separately with angular velocities w<. This critical current I c ,i has a remarkable 
non-monotonous dependence on the applied magnetic field which is due to a dynamically-induced 
structural transition with a rearrangement of vortices over the shells similar to the Coster-Kronig 
transition in hollow atoms. Thermally-activated externally-driven flux motion in a disk with pinning 
centers explains experimentally observed uii as a function of Io and T and the dynamically-induced 
melting transition. 

PACS numbers: 74.25.Qt, 74.25.Sv, 74.78.Na 



I. INTRODUCTION 

A mesoscopic Corbino disk is a unique system to 
study the dynamics of self-organized vortex matter in 
small-size superconductors. The interplay between the 
vortex- vortex and vortex-boundary interactions in meso- 
scopic superconductors leads to shape-induced giant vor- 
tex states pj , concentric shells of vortices and symmetry- 
induced vortex-antivortex "molecules" in mesoscopic 
squares and triangles 0- In the Corbino geometry Q 
an applied current creates a gradient in the current den- 
sity and thus the Lorentz force, i.e., introducing a shear 
driving force between the rings of vortices. This gives us 
the unique opportunity to study various dynamical ef- 
fects related to vortex motion, e.g., the transition from 
elastic to plastic motion, channeling [4|, vortex friction 
P|, etc. The dynamics of self-organized mortex matter 
in mesoscopic disks has many common features to, e.g., 
atomic matter, charged particles in Coulomb crystals, 
vortices in rotating Bose-Einstein condensates, magnetic 
colloids, synthetic nanocrystals, etc. 0, 0, Hi, or even 
large charged balls diffusing in macroscopic Wigner rings 
and can provide us with a deeper understanding of, 
e.g., the microscopic nature of friction, transport, mag- 
netic, optical and mechanical properties of various phys- 
ical and biological systems. 

In a Corbino disk, the applied current is injected at the 
center and removed at the perimeter (see Fig. 1) to induce 
a radial current density J that decays as \ j p along the 
radius As a result, vortices near the center of the 

disk experience a stronger Lorentz force Fl than those 
near the disk's edge. For small J, the local shear stress is 
small and the whole vortex pattern moves as a rigid body. 
Larger J result in a strong spatially inhomogeneous stress 
that breaks up the vortex solid and concentric annular 
regions move with different angular velocities. The volt- 
age profiles measured in experiments Q reflect different 
dynamical phases (elastic motion, shear-induced plastic 
slip) of vortex motion. The onset of plasticity in large 
Corbino disks was theoretically analyzed within a con- 



tinuous model in [Tl|. Within molecular dynamics (MD) 
simulations of interacting vortices at T = 0, the nuclc- 
ation and motion of dislocations in the vortex lattice was 
studied in |12j . 

Recently, using the Bitter decoration technique, the 
first direct observation of rings of vortices for mesoscopic 
superconducting Nb disks was reported |l|. For vor- 
ticities L = to 40, the circular symmetry led to the 
formation of concentric shells of vortices , sim ilar to elec- 
tron shells in atoms or in nano-clusters [l4|. The analy- 
sis of shell fil ling revealed "magic- number" configurations 
(MNC) pjj, [la, C3 corresponding to a commensurabil- 
ity between the shells which occurs when the numbers of 
vortices of each shell have a common divider . 

Here, we study the dynamics of vortex shells in meso- 
scopic Corbino disks. Our system has the added flexibil- 
ity that we have several experimentally accessible tuning 
parameters as, e.g., the driving shear force, the sepa- 
ration between the vortex rings (through the external 
magnetic field) and the commensurability between the 
vortex rings (through the relative number of vortices in 
each ring). In mesoscopic Corbino disks, we reveal a non- 
monotonous dependence of the critical current I c sep- 
arating two dynamical regimes, a "rigid body" or sep- 
arate shell rotation, on the magnetic field h, and the 
appearance of dynamical instabilities associated with a 
jump in I c (h). We show that this unusual behaviour is 
related to a "structural transition", i.e., an inter-shell 
vortex transition. For non-zero temperature, thermally- 
activated externally-driven flux motion is investigated, 
and we explain the observed two-step melting transition 
in Corbino disks |3|. 



II. THEORY AND SIMULATION 

We place a Corbino disk which has thickness d and 
radius R in a perpendicular external magnetic field Hq. 
The Corbino setup is shown in Fig. 1. An external cur- 
rent flows radially from the center to the edge of the disk 
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FIG. 1: The Corbino setup: the applied current is injected 
at the center and removed at the perimeter of the disk to 
induce a radial current density J (shown by dark blue/dark 
grey arrows) that decays as 1/p along the radius. As a re- 
sult, vortices near the center of the disk experience a stronger 
Lorentz force Fl (shown by yellow/light grey arrows) than 
those near the disk's edge. The vortices are shown by red- 
to-yellow/grey-to- light grey tubes and by red/grey spots on 
the surface. The direction of the external applied magnetic 
field Ho, which is perpendicular to the surface of the disk, is 
shown by dark yellow/grey arrow. 

and results in the inhomogeneous sheath current density 
J(p) = Iext/Zitp, which makes vortices closer to the cen- 
ter feel a stronger force compared to the ones near the 
edge. The Lorentz force (per unit length) acting on vor- 
tex i, <&oJ x z, resulting from the external current is: 
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where pi is the vortex position, = pi/R, and z 
is the unit vector along the magnetic field direction 
which is taken perpendicular to the disk. Here fo = 
= 47r n £, 2 Hi /R is the unit of force, I = 
/i A/ ext /$ , and A = X 2 /d. 

In a thin superconducting disk such that d < £ <C 
R <C A, the vortex-vortex interaction force /™ and the 
force of the vortex interaction with the shielding currents 
and with the edge ff can be modelled respectively by 
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where h = itR 2 ij i0 H /<S> = (H /2H c2 )(R/£) 2 , ^ is the 
position of the zth vortex, and L is the number of vor- 
tices, or the vorticity. Our numerical approach is based 
on the Langevin dynamics algorithm, where the time in- 
tergartion of the equations of motion is performed in the 
presence of a random thermal force. The overdamped 
equations of motion becomes: 



rjVi 



(4) 



Here ff — ff is the driving force (Eq Q), f£ p is the 
force due to vortex-pin interactions |21|. and ff is the 
thermal stochastic force, obeying the fluctuation dissipa- 
tion theorem 

< UiWfpjtf) >= 2 v 5 a(3 S lJ 6(t~t')k B T, (5) 

where Greek and italic indices refer to vector components 
and vortex labels, r\ is the viscosity. The ground state of 
the system is obtained by simulating field-cooled experi- 
ments 221. 



III. VORTEX DYNAMICS IN TWO-SHELL 
CORBINO DISK 

First, we consider the smallest mesoscopic Corbino 
disk which shows the main physics, and it has L = 19 
vortices which form the MNC (1,6,12), as shown in the 
left-hand inset of Fig. 3. At weak applied magnetic field 
h = 24 all the vortices are packed in an almost perfect 
triangular Abrikosov lattice. (Here we assume a perfect 
disk with no pinning, and temperature is set to zero after 
the "annealing" process.) We apply an increasing exter- 
nal current Jo, and we study the average angular velocity 
LOi of vortices in each shell. For small Jo, the shear pro- 
duced by the gradient of the Lorentz force is insufficient 
to break the vortex lattice. It produces only elastic de- 
formations in the lattice which rotates as a rigid body 
with angular velocity 

wrb = J/2tt < r 2 >, 

where < r 2 >= T 1 /L 1S the average square radius 

for the vortices in the disk. 

With increasing Iq, the Lorentz force gradient reaches 
a critical value at which the shells start to slide with re- 
spect to each other and rotate with different velocities 
u>i > u>2- We call this the "critical current" I Cl and for 
the case shown in Fig. 2(a) we found I c w 10.4. Now 
lui > ujrb > ^2, and the inner (outer) shell rotates faster 
(slower) than the rigid body. For higher applied magnetic 
field, e.g., h = 25, the critical current decreases, I c « 6.5 
(Fig. 2(b)). This is related to a deformation of the shells, 
which approach a circular shape, and therefore can more 
easily slide. At the same time the shells move closer 
to each other increasing the dynamical friction between 
them. At h = 26.25, the I c further decreases down to 
I c ~ 3, and the motion becomes unstable: for Jo > 3, 
the inner shell rotates with angular velocity u>2 > uii, al- 
though at some values of Jo in the region 3 < Jo < 6.8, 
u>2 drops down to 012 = wi (Fig. 2(c)). The appearance of 
this irregular motion, or "stochastization" , is related to 
the strong dynamical friction between the shells, which 
can even "lock in" the shells at some values of Jo . A fur- 
ther increase of the magnetic field, h = 26.5 (Fig. 2(d)), 
results in bistable motion of the shells: a second critical 
current J c i « 7.1 appears, and for currents I c < Iq < I c \ 
the system either rotates as a rigid body or the shells ro- 
tate separately. Finally, at even higher field, e.g., h — 30, 
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FIG. 2: The angular velocities u>i (orange/light grey solid 
circles) and ui2 (blue/dark grey open squares), measured in 
units of ljo = 'iivfioHc!; 2 /r]R 2 , of vortices in the first and 
second shells, respectively, as a function of applied current 
To, measured in units of /ioA7 ea; t/$o, for L = 19 and different 
magnetic fields h: (a) h — 24, a rigid-body rotation for Jo < 
I c ~ 10.4; the shells rotate separately for To > T c ; (b) h — 25, 
the critical current decreases, T c ~ 6.5; (c) h — 26.25, the T c 
further decreases, T c ~ 3, and the motion becomes unstable; 
(d) h — 26.5, bistable motion: a second critical current T c i ~ 
7.1 appears; (e) h = 30, the motion stabilizes, and first critical 
current T c disappears, and new critical current I c — I c ±. 
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FIG. 3: The critical current T c (in units of ^oh.I e xt/^o) versus 
h (in units of (T7o/2JT C 2)(-R/£) 2 ) near the structural transition 
between the states (1,6,12) (solid black circles) and (1,7,11) 
(solid red/grey triangles), insets. Empty symbols correspond 
to bistable vortex motion, with two critical currents, I c and 
Id- The jump in T c is related to an inter-shell vortex tran- 
sition from a higher to a lower orbit, similar to the Coster- 
Kronig transition in hollow atoms. 

the upper branch (corresponding to higher angular veloc- 
ity) of the function uj(Iq) for the inner shell (shown by 
orange/light grey dotted line in Fig. 2(d)) disappears. 
Therefore, the first critical current I c also disappears, 
and the motion stabilizes with the only critical current 
I c = Id (Fig- 2(e)). The non-monotonous change of the 
critical current I c with increasing applied magnetic field 
h is summarized in Fig. 3. First, I c decreases, then insta- 
bilities develop in the system which in general are indica- 
tive of, and precede, a phase transition. The transition 
to a second critical field occurs through a bistable state 
characterized by two critical currents in the system. This 
unusual behaviour can be understood by analyzing the 
critical current together with the structure of the vortex 
shells in the region of this "phase transition" . The vortex 
patterns for magnetic field before (low fields) and after 
(high fields) the transition are presented in the left-hand 
side and in the right-hand side insets of Fig. 3. The jump 
in the critical current is correlated to a "structural tran- 
sition" in the system where the distribution of vortices 
over the shells is altered, i.e., there is an inter-shell vor- 
tex transition from a higher orbit to a lower one, which is 
similar to the Coster-Kronig transition in hollow atoms 
(see, e.g., |23j). This structural transition is accompanied 
by a local re-distribution of the flux inside the disk, or 
by the appearance of flux jump instabilities, which might 
have an additional triggering mechanism caused by the 
viscosity. 

We found that this jump in I c is also observed for other 
vortex configurations and is generic for vortex structures 
confined in Corbino disk. For instance, in a disk with 
L = 22 it is related to the transition between the states 
(1,8,13) and (2,8,12) at h = 33.5. Note that a similar 
behavior (a jump in the mobility) is also found in the 
Frenkel-Kontorova model for the locked-to-sliding tran- 
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FIG. 4: The angular velocities u>i (in units of u>q = 
47r (ioH 2 £ 2 1 'r/R 2 ) versus Iq (in units of fioJ^-Iext/^o) in first 
(orange/light grey open circles), second (black solid line) 
and third (blue/dark grey open squares) shells, for a three- 
shell vortex configuration with L — 37. The inset shows 
a triangular-lattice ground-state configuration (1,6,12,18) for 
h = 37 at Io = 0. The critical currents are: I c \2 ~ 9.1, and 
/ c23 ~ 9.5. 

sition for chains moving on commensurate potentials |24| . 

IV. THERMALLY-ACTIVATED 
EXTERNALLY-DRIVEN DYNAMICS IN A 
CORBINO DISK 

Consider now a Corbino disk containing a larger num- 
ber of vortex shells, e.g., L — 37 where vortices form 
three shells. Note that the chosen vorticity allows the 
triangular-lattice MNC (1,6,12,18) shown in the inset of 
Fig. 4. As in the above case of two shells, for small Io 
the system displays a rigid-body rotation. However, the 
distance between the inner and the outer shells is larger, 
and, thus, the Lorentz force gradient results in a stronger 
elastic deformation of the vortex configuration as com- 
pared to the two-shells case. The shear stress is stronger 
close to the center of the disk, and as a result, the in- 
ner shell splits off first at J c i2, with increasing Jo, while 
second and third shells still rotate with the same angular 
velocity up to some critical current value I C 23 > I c i2 (see 
Fig. 4). 

The three-shell system displays a number of inter-shell 
vortex transitions when changing the applied magnetic 
field. Similar to the case of two shells, a structural "phase 
transition" occurs where vortices from a higher shell tran- 
sit to a lower shell when the magnetic field increases. 
This leads to non-monotonic dependences of the critical 
currents I c \2 and I C 23 on h and to instabilities in the 
functions u>i{Io) at different h. 

Now we will vary the temperature T and study the 
thermally-activated externally-driven dynamics of vor- 
tices in the shells. For this purpose, we introduce pinning 
centers randomly distributed over the disk (for brevety 
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FIG. 5: The angular velocities of vortices LUi (in units of too = 
A-khoH 2 ^ 2 /nR 2 ) versus T (in units of T = Air^Hl^d/kBs) 
in first, second and third shells (using same symbols as in 
Fig. 4), for L = 37, I = 0.4, f p /f Q = 2.5. The onset of 
motion occurs in the form of rigid-body rotation. At T = T* lt 
the inner shell splits off, the second shell unlocks at T = T* 2 i , 
and shells rotate independently ( angular melting) . The radial 
melting, i.e., dissolving of the shells, occurs at T — Tm,t- The 
inset: the experimental T-dependence of LUi obtained from the 
resistivity measurements at three different distances from the 
center of a Corbino disk 

we assume dense narrow pinning centers with maximum 
pinning force / pU ) and apply a very low current Iq that 
results in a Lorentz force which is less than the pinning 
force f p o for a vortex in any shell. At low T all the 
vortices are pinned and the shells do not move. With in- 
creasing T, the random-noise force is added to the (weak) 
Lorentz force. Intuitively one expects that vortices in the 
inner shell, where the Lorentz force is maximum, will de- 
pin first and start to rotate while vortices in the other 
shells are still pinned. However, the vortex-vortex inter- 
action prevents this scenario, and the inner-shell vortices, 
if unpinned, become "trapped" by the vortex lattice itself 
that leads to an elastic deformation. The vortex-vortex 
interaction locks in the vortex configuration, and when 
for further increasing T, the onset of the motion occurs in 
the form of a rigid-body rotation (Fig. 5). At some tem- 
perature T* 1; the inner shell splits off and starts moving 
with a higher angular velocity, while higher shells stay 
locked up to a "second critical" temperature T* 2 . The 
transition from a rigid body to a separate rotation of 
shells is called angular melting (i.e., the "vortex solid" 
to "vortex shells" transition). This is a multi-step pro- 
cess which starts when first shell splits off (i.e., at T* x 
in our three-shell model) and it finishes when all the 
shells rotate independently (T > T* 2 ), but they keep 
their identity and contain a constant number of vortices 
with well-defined average radius {rf). The radial melting 
is associated with a dissolving of the shells (i.e., the "vor- 
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tex shells" to "vortex liquid" transition) and occurs at a 
higher temperature, Tm, t (Fig- 5). 

This three-shell model explains qualitatively the ex- 
perimentally observed thermally-induced "solid-liquid" 
melting transition in a Corbino disk 3] with three probes 
(see inset to Fig. 4). Also, the behaviour of Wi(Io) (Fig. 4) 
is in agreement with the experimental results. In the 
experiment, the superconducting disk contained a large 
number of vortices, and the measurements were done for 
rings rather than for shells. To model larger systems, we 
used a disk with, e.g., 123 vortices forming six shells, and 
we calculated average velocities in three rings (each con- 
taining two shells). We found that the results for larger 
disks are similar to the above three-shell system, but the 
radial melting temperature Tm,t is lower, in agreement 
with the experiment. 

V. CONCLUSIONS 

We predict an unusual non-monotonous behavior of 
the critical current for unlocking of the vortex rings with 



magnetic held which originates from a "structural transi- 
tion" where a vortex jumps from a higher shell to a lower 
shell similarly to the Coster-Kronig transition in hollow 
atoms. Pushing the experiments of Ref. into the meso- 
scopic regime will allow one to detect experimentally the 
predicted behaviour. The vortex motion in the presence 
of pinning centers reveals the onset of a rigid body ro- 
tation, due to thermally-activated externally-driven flux 
motion. With increasing T, the inner shell splits off first 
and subsequently all the shells start moving separately. 
The present numerical study explains the experiments of 
Ref. |3J, in which vortex velocities in adjacent layers in a 
Corbino disk were studied as a function of / and T. 



ACKNOWLEDGMENT 

We thank B. Janko, F. Nori, and U. Welp for useful 
discussions and hospitality. This work was supported by 
the Flemish Science Foundation, and the Belgian Science 
Policy. V.R.M. acknowledges a partial support through 
POD. 



[1] V.A. Schweigert, F.M. Peeters, and P. Singha Deo, Phys. 
Rev. Lett. 81, 2783 (1998); A. Kanda, B.J. Baelus, 

F. M. Peeters, K. Kadowaki, and Y. Ootuka, Phys. Rev. 
Lett. 93, 257002 (2004). 

[2] L.F. Chibotaru, A. Ceulemans, V. Bruyndoncx and V. V. 
Moshchalkov, Nature (London) 408, 833 (2000); Phys. 
Rev. Lett. 86, 1323 (2001); V.R. Misko, V. M. Fomin, 
J. T. Devreese, V. V. Moshchalkov, Phys. Rev. Lett. 90, 
147003 (2003). 

[3] G.W. Crabtree, D. Lopez, W.K. Kwok, H. Safar, and 

L.M. Paulius, J. Low Temp. Phys. 117, 1313 (1999); 

D. Lopez, W.K. Kwok, H. Safar, R.J. Olsson, A.M. Pe- 

trean, L.M. Paulius, and G.W. Crabtree, Phys. Rev. Lett. 

82, 1277 (1999). 
[4] R. Besseling, P.H. Kes, T. Drose, and V.M. Vinokur, New 

J. Phys. 7, 71 (2005). 
[5] A. Maeda, Y. Inoue, H. Kitano, S. Savel'ev, S. Okayasu, 

I. Tsukada, and F. Nori, Phys. Rev. Lett. 94, 077001 

(2005). 

[6] T.B. Mitchell, J.J. Bollinger, W.M. Itano, and 
D.H.E. Dubin, Phys. Rev. Lett. 87, 183001 (2001). 

[7] A. Pertsinidis and X.Z. Ling, Nature 413, 47 (2001). 

[8] C. Murray, C. Kagan, and M. Bawendi, Ann. Rev. Mat. 
Sci. 30, 545 (2000). 

[9] C. Coupier, M. Saint Jean, and C. Guthmann, Phys. Rev. 
E 73, 031112 (2006). 
[10] Y. Paltiel, E. Zeldov, Y. Myasoedov, M.L. Rappaport, 

G. Yung, S. Bhattacharya, M.J. Higgins, Z.L. Xiao, 
EY. Andrei, P.L. Gammel, and D.J Bishop, Phys. Rev. 
Lett. 85, 3712 (2000). 

[11] P. Benetatos and M.C. Marchetti, Phys. Rev. B 65, 
134517 (2002). 



[12] M.-C. Miguel and S. Zapperi, Nature Materials 2, 477 
(2003). 

[13] I.V. Grigorieva, W. EscofHer, J. Richardson, LY. Vin- 
nikov, S. Dubonos, and V. Oboznov, Phys. Rev. Lett. 
96, 077005 (2006). 

[14] Y.N. Ovchinnikov and V.Z. Kresin, Eur. Phys. J. B 45, 
5 (2005); 47, 333 (2005). 

[15] V.A. Schweigert and F.M. Peeters, Phys. Rev. B 51, 770 

(1995) . 

[16] V.M. Bedanov and F.M. Peeters, Phys. Rev. B 49, 2667 
(1994). 

[17] L.J. Campbell and R.M. Ziff, Phys. Rev. B 20, 1886 
(1979). 

[18] A.I. Buzdin and J. P. Brison, Phys. Lett. A 196, 267 
(1994). 

[19] L.R.E. Cabral, B.J. Baelus, and F.M. Peeters, Phys. Rev. 

B 70, 144523 (2004). 
[20] B.J. Baelus, L.R.E. Cabral, and F.M. Peeters, Phys. Rev. 

B 69, 064506 (2004). 
[21] C. Reichhardt, C.J. Olson, and F. Nori, Phys. Rev. Lett. 

78, 2648 (1997); Phys. Rev. B 57, 7937 (1998); Phys. 

Rev. B 58, 6534 (1998). 
[22] K. Harada, O. Kamimura, H. Kasai, T. Matsuda, 

A. Tonomura, and V.V. Moshchalkov, Science 274, 1167 

(1996) . 

[23] J. Limburg, J. Das, S. Schippers, R. Hoekstra, and 
R. Morgenstern, Phys. Rev. Lett. 73, 786 (1994). 

[24] O.M. Braun and Yu.M. Kivshar, The Frenkel-Kontorova 
model: concepts, methods and applications (Springer, 
Berlin, 2004). 



